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Abstract 

Within the Matsubara formalism we derive expressions for the pion selfenergy and 
^ \ the two-pion propagator in a hot pion gas. These quantities are used to selfconsis- 

tently calculate the in-medium tttt amplitude beyond the quasiparticle approxima- 
tion (QPA). The results are shown to differ significantly from QPA-based calcula- 
tions. We also examine the impact of chiral constraints on the tttt interaction in a 
$— i ' chirally improved version of the Jiilich tttt model. 
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1 Introduction 



It is anticipated that future experiments at the Brookhaven Relativistic Heavy Ion Collider 
(RHIC) will create a highly excited, but nearly baryon-free central zone. At later stages 
of the collision this zone will contain a dense meson gas, primarily consisting of pions 
with temperatures of the order of the pion mass. The equation of state and the transport 
properties of the pion gas are therefore of great current interest Jl|, |2|. One expects 
that the nir interaction will be modified due to statistical correlations (Bose occupation 
factors) as well as dynamical (selfenergy) effects |3], [|, ^, [7]. In this context we have 
previously presented a Brueckner-type calculation for a thermal pion gas treating both 
effects selfconsistently ||. To facilitate the numerical calculations we have employed 
the quasiparticle approximation (QPA) for the single-pion propagator. We have noticed, 
however, that this approximation becomes unreliable for temperatures above T w 150 
MeV, i.e. in the vicinity of T c , the chiral restoration temperature. In addition, recent 
studies of pion correlations in cold nuclear matter |J have shown that, for a thorough 
description of the subthreshold region, it is mandatory to work with a tttt model which 
obeys constraints imposed by chiral symmetry. Since chiral symmetry predicts properties 
of the 7T7T interaction in the soft pion limit (m n — »0), there is an intimate relation between 
the proper off-shell treatment of in-medium properties and the implementation of chiral 
constraints. 

The aim of the present paper is twofold: first we investigate effects of proper off-shell 
calculations by comparing them to the QPA results of ref. ||. This will be discussed in 
sect. 2. To incorporate constraints from chiral symmetry we then present in sect. 3 an 



improved version of the Tin interaction of Lohse et al. |L0|, [Tl[] which is subsequently used 
for the in-medium calculations. A short summary and concluding remarks are given in 
sect. 4. 



2 



2 Formalism and Results from the Jiilich nn Model 



In the 7T7T meson exchange model of Lohse et al. (T^, |TT| the vacuum tttt T-matrix in 



a 



given partial- wave/isospin channel, JI, is obtained from a Lippmann-Schwinger equation, 



where the Blankenbecler Sugar reduction scheme |12] of the 4-dimensional Bethe Salpeter 
equation is adopted: 

T^{E,q x ,q 2 ) = V^{E,q x ,q 2 ) 

oo 

+ J dk k 2 Au 2 k V£(E, qi, k) G n7T (E, k) T£(E, k, q 2 ) . (1) 
o 

Here E denotes the starting energy of the pion pair and q\, q 2 are the relative CMS 3- 
momenta of one pion in the in- and outgoing state, respectively. The interaction kernel 
is calculated from an effective meson Lagrangian including p-, / (1400)- and /2(1270)- 
s-channel exchanges, a p-meson exchange in the t-channel as well as coupling to the KK 
channel (not written explicitely in eq. (1)). The vacuum BbS two-pion propagator in the 
CMS reads 

G°^(E,k) = j 1 -^- D^k) D n (E-k ,k) 
11 

~ k E 2 - 4ul + irj ' (2) 
where D^^ko, k) denotes the single-pion propagator and uik = m 2 + k 2 . 
At moderate gas densities it is sufficient to treat medium modifications at the one- and 
two-body level rendering a Brueckner-type description. Leaving the pseudopotentials 
temperature independent, which should be a good approximation for temperatures well 
below the masses of the exchanged mesons, the temperature dependence of the T-matrix 
is solely induced by the thermal two-pion propagator, CrI ' . In terms of the thermal 



single-pion propagator, D^, it can be evaluated by means of a Matsubara sum |T3|, [TjJ : 

C^(0 A , h, k 2 )=T ]T D v (z v , h) D n (n x - z u , k 2 ) , (3) 
with discrete frequencies 

z v = Q\ = u v = 1-kvT , u, X integer. (4) 



Using standard techniques of analytical continuation one then obtains 

GUE + , h, k 2 ) = \( ^^I m D n (u + , h) ImD w (u;' + , k 2 ) (1 + /» + />')) 
L J ix ix 

X f 7 77- + " 77-] (5) 

where u± = uo etc. and f n {uj) = (exp[o;/T] — l) -1 . For the T-matrix evaluation it 
will be sufficient to consider for E>0 and in the two-pion CMS, i.e. in back-to-back 
kinematics with &2=— k±. In this case the imaginary part of eq. (5) receives 2 contributions: 



duo , 



ImG^ A (E + , k) = — f — [1 + r(u) + r(E - u)]ImD 7r (u+, k)ImD w {E - k), 

J 7X 


oo 

ImG^ B (E+, k) = - I —[r{LO — E) — r(u)]ImD n (u+, k)ImD n (u + -E,k) . (6) 

J 7X 



In the zero-temperature limit contribution A reproduces the imaginary part of eq. (2), 
whereas B vanishes. The latter arises solely from thermal excitations. The real part of 
G^ n is calculated via a dispersion integral B: 

o 

In a second step we evaluate the pion selfenergy, £„., which enters the single-pion propa- 
gator 

At the two-body level H n is expressed in terms of the txtx invariant forward scattering 
amplitude as 

=- / T Y\M^{z v + z v ,,k,p)D 1 ,{z v >,p) , (9) 

J ( 27r ) 

which involves a Matsubara sum over z u /. Using the spectral representations of M n7T and 
D n one arrives at [|15[ 



H w {u + ,k) = / -—r ImM n7r (E +J k,p)ImD 7T (uj + ,p) ■ — — . (10) 

J ylix) 6 ix ix uo + + to — E 
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The imaginary part of £„- is given by the 5- function contribution of l/(u + + u' — E'). The 
positive-energy part of the u' integration gives rise to the standard contribution commonly 
used in the literature: 

(11) 

Note that the Matsubara formalism yields the 'correction' f n (u> + a/) in the occupation 
factor as was found by Chanfray et al. |7| in the real-time formalism. The negative-energy 
part of the u' integration can be rewritten by exploiting the symmetry properties of the 
1-pion propagator, the M^-amplitude and the Bosefactor resulting in 

J —ImMm(\u-cJ\ + ,k,i?)ImD n (u' + ,p) 

[sgn(u-u')r(u')-r(\u-u'\)) . (12) 

with sgn(u; — uf) = ±1. In analogy to the thermal two-pion propagator (eq. (6)) the con- 
tribution B accounts for thermal excitations of the system. E.g. in lowest order, JmS^ 
is only non- vanishing for uo > 3m,r. The M^-amplitude is related to the T^-matrix 
after transforming it into the CMS and identifying the starting energy as E 2 = s = 
(lu±lu') 2 -(p + k) 2 : 

M n7T (uj ± u', k, p) = M^(v^, q, q) = (2tt) 3 4lu 2 T nn (E, q, q) (13) 

with q = q(u,u', k,p, cos 9) being the CMS momentum of both in- and outgoing pions. 
After a variable transformation cos© — > E, where = l(k,p), the final expressions read 

ImE 7r (oj + ,k) = JmEj((i) + , k) + JmSf k) 

IrnSi( U+ ,k) = -\]^J d ^ImD^,p)[r{u')-r{u + u')]I A 

o ^ ' o 

ImS%(u, + ,k) = -lJ^J^/ m ^(a;;,p)[ S gn(a;-a;0/V) + r(l^-^l)]^ 



(14) 
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with 



„A/B 
rn.a :r 



I A/B = J dE E M nn (E,q,q) , (15) 
where the upper and lower bounds, 



A/B 
min 



Emax/min = (^ W) - k -p ± 2kp 

J max /min 



ELxImin = (" - u/) 2 -k 2 -p 2 ± 2k P , (16) 



correspond to cos = =pl, respectively. For the numerical evaluation of ImM^ we include 
all partial waves up to J = 2 and all isospins. Eq. (10) evades a 'direct' calculation 
of i?eS 7r (co>, k) in terms of ReM vlT due to the .E'-dependence of the occupation factor. 
Therefore we again make use of a dispersion relation including a subtraction at u = 0: 

Re^k) = -Vj Q — -p. (17) 

The subtraction is dictated by the fact that ReH^uo^k) — > for u; — > 0. To check 
the subtraction procedure numerically we have performed calculations by approximating 
f n (E') ~ f n (w + w') in eq. (10), which allows to reinsert the spectral representation 
of M^. ReYiK can then be calculated directly in terms of ReM n7T (in exact analogy to 
calculating ImH n in terms of ImM^). The results of this test calculation coincide within 
a few percent with the results obtained from eq. (17). Eqs. (1), (6)-(8) and (14), (17) 
form the selfconsistent set of equations which is solved by numerical iteration. 
The results for the imaginary part of the 7T7t M- amplitude are displayed in Fig. 1. With 
increasing temperature we find in the cx-channel ( J/=00) a considerable enhancement over 
the vacuum curve for CMS energies below E^450 MeV. Apart from the highest tempera- 
ture (T=200 MeV) there is, however, no significant strength below Im^. For T=200 MeV 
some accumulation of subthreshold strength is found which is almost entirely due to ther- 
mal excitations (the 'B '-contribution in eq. (6)). For comparison, results within the QPA 
are also given (in contrast to ref. || the 'polestrength' Zk = (1 — dRell^oj, k) / duj 2 ^)^ 1 
is not included since it leads to a violation of the normalization condition of the 1-pion 
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spectral function). For the p-channel {J I =11) both the full and the QPA results show a 
broadening of the resonance, increasing with temperature due to the thermal motion of 
the gas particles. The QPA results clearly overestimate the thermal broadening and the 
upward peakshift of the p resonance. 

A comment concerning the QPA calculations of refs. || |H| is in order. When expanding 
around the quasiparticle pole, in next-to-leading order a complex polestrength z k arises. 
This results in a strong enhancement of ImM n7T in the vicinity of resonance peaks. Such 
a strong enhancement is not confirmed by the full calculations presented here. At in- 
termediate energies the full calculations do show an enhancement over the QPA results, 
however. 



3 7T7T Interaction with Chiral Constraints 

For the subthreshold irn amplitude in cold nuclear matter the implementation of properties 
dictated by chiral symmetry has turned out to be crucial 0. To see the impact of such 
constraints on the interacting pion gas we have improved the Julich model employed in 
the previous section in two respects (for a detailed discussion see refs. fl7|. |§ ) : 



- we introduce rnr contact interactions as obtained from the gauged nonlinear a model 
of Weinberg |18| . At the tree level these contributions ensure the correct values of 



the 7T7T amplitude at both the soft point and the Adler point. When supplemented 
with phenomenological form factors, F(q) = (2A 2 — 4m 2 ) 2 /(2A 2 + 4g 2 ) 2 , they are 
suitable for iteration in the T-matrix equation (1); 

instead of using the on-energy-shell prescription implied by the BbS formalism we 
here employ an on-mass-shell prescription for the off-shell values of the pseudopo- 
tentials V^J . In this way we are able to impose the correct chiral limit of the s-wave 
scattering lenghts, i.e. a 01 — > for m n — > (the model of Lohse et al. |10| does not 
have this property). 
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With similar cutoff parameters as in |T(J the results of refitting the vacuum tttt scattering 
data are quite satisfactory (see Fig. 2). The additional cutoff parameters in the contact 
terms are chosen to be slightly below 1 GeV. For the s-wave scattering lengths we obtain 
a 00 = 0.21m" 1 , a 02 = — 0.028m" 1 , in reasonable agreement with a recent experimental 
analysis [[19] . 



With this 'chirally improved' model we have repeated the selfconsistent in-medium cal- 
culations described in the previous section. The results are displayed in Fig. 3. We again 
observe a reduction of the peak values in both a- and p- channel of the tttt scattering 
amplitude. The effects are smaller than in the original Jiilich model. There remains 
some enhancement over the vacuum curve in the low-energy region and there is very little 
strength below threshold. This can be attributed to the subthreshold repulsion induced 
by the contact interactions. 

4 Summary 



Within the finite temperature Green's function method |13|, [15| we have derived expres- 
sions for the pion selfenergy and the uncorrelated 2tt propagator in a hot pion gas in 
thermal as well as chemical equilibrium (p^ = 0). Starting from the vacuum titt model of 



Lohse et al. |T0] these have been used to calculate medium modifications at the one- and 



two-body level selfconsistently, i.e. the in-medium scattering amplitude was evaluated 
in terms of the pion selfenergy and vice versa. It has been shown that the full off-shell 
calculations differ from the QPA in the following respects: in the intermediate-energy 
range (E^300-800 MeV) the imaginary part of the titc amplitude shows less suppression 
compared to the vacuum curve; thermal excitation contributions (not present in the QPA) 
lead to some subthreshold strength in the a channel at temperatures around 200 MeV; 
in particular, the QPA calculations overestimate the thermal broadening and the upward 
peakshift of the p resonance. 

Motivated by a recent analysis of titt correlations in cold nuclear matter 0] , we have also 
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studied a 'chirally improved' version of the Jiilich tttt model. Applying this model to the 
hot pion gas it is found that the subthreshold strength almost entirely disappears due 
to the additional repulsion induced by the chiral constraints. Also the enhancement just 
above the two-pion threshold is much less pronounced. The in-medium features of the 
tttt amplitude for higher energies (reduction of the peak values in both a and p channel, 
essentially no peakshift of the p resonance) are in good agreement for both models. After 
all, the deviations from the nir interaction in free space are rather small. 
In conclusion, we wish to stress that proper off-shell calculations beyond the QPA and the 
inclusion of constraints from chiral symmetry are necessary to obtain reliable results for 
the pion selfenergy and the in-medium tttt amplitudes of a hot pion gas. The latter could 
serve as an input for transport models simulating the time evolution of ultrarelativistic 



heavy- ion collisons pOfl near freeze-out. 



Acknowledgement 



We thank G.G. Bunatian, G. Chanfray, J. Durso, B. Friman and J. Speth for useful 
discussions. One of us (R.R.) acknowledges financial support by Deutscher Akademischer 
Austauschdienst (DAAD). This work was supported in part by a grant from the National 
Science Foundation, NSF-PHY-89-21025. 



9 



References 

[1] G.G. Bunatian and B. Kampfer, preprint FZR 93-28. 

[2] S. Gavin, Nucl. Phys. B351 (1991) 561. 

[3] E.V. Shuryak, Nucl. Phys. A533 (1991) 761. 

[4] Z. Aouissat, G. Chanfray, P. Schuck and G. Welke, Zeit. Phys. A340 (1991) 347. 
[5] A. Schenk, Nucl. Phys. B363 (1991) 97. 

[6] H.W. Barz, G.F. Bertsch, P. Danielewicz and H. Schulz, Phys. Lett. B275 (1992) 
19. 

[7] G. Chanfray and D. Davesne, to be published. 

[8] R. Rapp and J. Wambach, Phys. Lett B315 (1993) 220. 

[9] Z. Aouissat, R. Rapp, G. Chanfray, P. Schuck and J. Wambach, Nucl. Phys. A581 
(1995) 471. 

[10] D. Lohse, J.W. Durso, K. Holinde and J. Speth, Phys. Lett. B234 (1989) 235; 
Nucl. Phys. A516 (1990) 513. 

[11] B. C. Pearce, K. Holinde and J. Speth, Nucl. Phys. A541 (1992) 663. 

[12] R. Blankenbecler and R. Sugar, Phys. Rev. 142 (1966) 1051. 

[13] L.P. Kadanoff and G.A. Baym, 'Quantum Statistical Mechanics', Benjamin, New 
York 1962. 

[14] L.D. Fetter and J.D. Walecka, 'Quantum Theory of Many Particle Systems', Mc 
Graw Hill, New York 1971. 

[15] M. Schmidt, G. Ropke and H. Schulz, Ann. Phys. 202 (1990) 57. 



10 



[16] R. Rapp and J. Wambach, Nucl. Phys. A573 (1994) 626. 

[17] J. Durso, R. Rapp and J. Wambach, to be published. 

[18] S. Weinberg, Phys. Rev. Lett. 17 (1966) 616 

[19] H. Burkhardt and J. Lowe, Phys. Rev. Lett. 67 (1991) 2622. 

[20] G.M. Welke and G.F. Bertsch, Phys. Rev. C45 (1992) 1403; H.W. Barz, P. 
Danielewicz, H. Schulz and G.M. Welke, Phys. Lett. B287 (1992) 40. 



11 



FIGURE CAPTIONS 



Fig. 1: The imaginary part of the irir M-amplitude in a (7 7=00) and p (77=11) channel 
in a hot pion gas at p w =0. The long-dashed, short-dashed and dotted lines are the 
selfconsistent results of the full off-shell calculations at temperatures T=150,175 and 
200 MeV, respectively. The dashed-dotted line represents the selfconsistent QPA 
calculation at T=200 MeV. The full line corresponds to the vacuum amplitude 
within the non-chirally symmetric Jiilich model. 

Fig. 2: A fit to the vacuum mi scattering phase shift data obtained within the chirally 
improved Jiilich model. 

Fig. 3: The imaginary part of the nn M-amplitude in a and p channel in a hot pion gas 
at /i n —0. All curves are calculated by means of the chirally improved Jiilich model 
(full line: free space; temperature identification as in fig. 1). The in-medium results 
are from selfconsistent, full off-shell calculations as described in sect. 2. 
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